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ABSTRACT 


The purpose of the investigation is to find how the natural 
frequency of transverse vibration is affected when the center of gravity 
and center of stiffmess are not coincident and to determine if it is 
important in calculating the natural frequency of the fundamental mode 
of vibration of ships. 


Two methods, one an approximate analytic method, and the other 
a lumped mass method, are used to determine the strength of the coupling 
of these modes of vibration. The approximate analytic method involves 
the solution of differential equations from an assumed mode shape and 
then the use of Stodola's equation to compute the change in natural 
frequency. The lumped mass method is a numerical procedure requiring a 
considerable amount of computation. Both methods are in good agreement 
when the ratio of natural transverse to longitudinal frequency is less 
than one-half. 


The approximate analytic method shows that the change in 
transverse frequency of vibration is a function of the product of two “ 
factors. One factor is directly proportional to the area times the distance 
of the center of gravity from the center of stiffness squared and inversely 
proportional to the area moment of inertia. The other factor is a com- 
plicated expression, but which is dependent only on the ratio of trans- 
verse to longitudinal frequency as computed by simpler theory in beams, 
neglecting the offset of the center of gravity from the center of stiffness. 
These expressions have been evaluated and plotted against frequency ratio. 
A simple method of computing the change in transverse frequency of vibration 
is provided. 


When the important parameters for calculating the change in 
frequency of vibration are calculated for a merchant ship and submarine, 
it is seen that the coupling of the fundamental modes of vibration in 
these structures is insignificant. 


In the higher modes of vibration, the ratio of transverse to 
longitudinal frequencies approaches unity. It is recommended that the 
coupling in these higher modes be investigated. 


Thesis Supervisor: Frank M. Lewis 
Title: Professor of Marine Engineering 
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NOTATION 


area of a cross section which will support longitudinal stress 
arbitrary constant 
modulus of elasticity 


ratio of transverse to longitudinal frequency computed disregarding 
the offset of the center of gravity from the center of stiffness 


force 

gravity 

area moment of inertia of a cross section 
spring constant 

length of a beam 

distance between masses in the lumped mass method 
moment 

mass 

stress 

longitudinal displacement 

shear 

frequency 


longitudinal frequency as computed disregarding the offset of the center 
of gravity from the center of stiffness. 


transverse frequency as computed disregarding the offset of the 
center of gravity from the center of stiffmess. 


transverse frequency with coupled longitudinal motion 
vertical displacement 

distance of the center of gravity from the center of stiffmess 
Slope 


a number determined by the transverse mode of vibration. It equals 
4.73 for the fundamental mode. 


weight per unit length 
-Vie 





Ris 


INTRODUCTION 


The calculation of the natural frequencies of the longitudinal and 
transverse modes of vibration of any structure have long been considered 
Mraebendent. \>) In a uniform bar, of course, it is true, but for example, 
in a ship's hull where the center of mass and center of stiffness are not 
coincident, it is not necessarily true. A bending of the ship's structure 
causes some motion in the longitudinal direction of the center of gravity. 
This, in turn, will cause some change in the natural frequency of transverse 
vibration. It is the purpose of this thesis to investigate the strength of 
the coupling of longitudinal and transverse modes of vibration when the 
centers of gravity and stiffness are displaced and determine if it has any 


significance in the calculation of natural frequencies of ships. 


5. 
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o—_ <= - = —“2 
—_ <j 66 
- =< — 
a _ 
7 6) a> ~ 
a = -= _> © £2 =| 
a = 





== 
PROCEDURE 


Two distinct procedures are used to find the coupling of longitudinal and 
flexural modes of vibration; one method is an approximate analytic method 
and the other is a lumped mass method. The steps to each method are outlined 
in subsections A and B below; the details are in the appendices. Through 
a numerical example the results of the two procedures are then compared 
for agreement. The entire analysis is then interpreted in the light of 
typical values of significant parameters for a submarine and a merchant ship. 

The investigation is restricted to the fundamental mode of 
vibration. 

In neither the approximate analytic method nor the lumped mass 
method is the deflection of the system due to shear and rotary inertia 
considered. The transverse motion is due to bending only. 

A. Approximate Analytic Method. 

The general procedure in this method is to find a correction to an 
assumed mode shape from the differential equations of a system in which 
bending causes longitudinal motion and longitudinal forces cause bending. 
The change in mode shape is then used as a measure of the change in natural 
frequency of transverse vibration. The steps in the procedure which are 
carried out in Appendix I are: 

1. Determine the differential equations of motion of a system in 

which the center of gravity and center of stiffness are not 
coincident. The first equation is a longitudinal force balance 


equation and the second one is a bending moment equation. 


2 2 
- (ah) Om (u + oy y) + AE dau = 0 (1) 
gz at Ox ax 





tof) 3° + EI ay + za 2 y =O (2) 
dt Ox ox 
(aA) = weight per unit length 
A = area of cross section 
g = gravity 
u = longitudinal displacement 
y = vertical displacement 
E = Young's modulus 
y = distance of center of gravity from the center of 
stiffness 
a = area moment of inertia of a cross section 
2) Assume a transverse mode shape of a free-free beam in space. 


3) 


4) 


>) 


y = ¢| cos x + (eA 5) Cosh X x | (3) 


Q 
ll 


arbitrary constant 


AL = @ number determined by the mode of vibration. It is 
equal to 4.73 for the fundamental mode of vibration. 


L length of the structure. 


From the differential equation for longitudinal force balance, 
find the longitudinal motion of the center of stiffness pro- 
duced by the mode shape. 

This longitudinal motion causes longitudinal forces which 
produce moments about the center of gravity and hence additional 
bending of the system. Determine the additional deflection 
from these forces using a bending moment equation. 

With the results of Step 4 recycle Steps 3 and 4. Hence, a 
first axa a@ second approximation to the additional deflection 


can be obtained. 
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6) Use a method similar to that of Stodola to find the trans- 
verse frequency with coupled longitudinal motion. The ratio 
of the square of the transverse frequency with longitudinal 
motion to that without is proportional to the ratio of total 
deflection without longitudinal motion to the total de- 


flection of the beam with longitudinal motion. 


B. Lumped Mass Method. 
The lumped mass procedure is the Prohl-Myklestad sequence method 


for calculating the natural frequency of beams of variable cross section. 
This calculation is modified to include the effect of the displacement of 
the center of gravity from the center of stiffness. 

The calculation is based on the fundamental equations of bending 


and strain in beams: 


= 

M da 

ine 5 and 
dx 

F - WE), 


The beam is divided into sections. At each section an inertia loading is 
added: 


in shear equal to 


2 
mW y 


in bending equal to 
mw” ¥ (u - ¥ 6) 

and in longitudinal force equal to 
m wo (u-y 0) 


m = mass 

w = frequency 
8 = slope 

£ = length 





=5= 


The procedure is to calculate the longitudinal force and displacement, trans- 
verse deflection, slope, bending moment and shear at each section along the 
beam in proceeding from the end of the beam to the mid-point. Various fre- 
quencies are tried until the boundary conditions are satisfied. For the 
fundamental mode, the slope, shear, and longitudinal motion at the mid-point 


are zero. 
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RESULTS 


A. Approximate Analytic Method 


the cquation for determining the frequency of transverse vibration 


by the approximate analytic method is: 


2 2 re) 
“= %p POPES (14) 
t ot Y_ + Y + Ys 
where 
7 2 f2ek PF Bet 
ot. . (oA) (5) 
: Cos (AL/2 
= el 2 Cos (AL/2) - 1 = aoaae Gs Cie (6) 


L/2 
aS ae Cos VD, 
Y i y c Cos RL / 2 pas oe 


, {Cos ae - 1) . il (Cosh AL/2 - 1) 


(Seis) comt/s (ebay et (7) 
é aps" ,' k D,* w 
Y,= (2) 7 Cos AL/2 i ) 
2 GF L y ¢ Oo ie aa oes 


(1 as 


L/2 " 
his D, ae aL Ae | a mA Sin 
ii L/2 - @ | “7 2 
> 2, 1 = Dd, Xd ca DS (Cos D, ) 


fi i (Cos aet2 LINN i aan 
it ih 
.- (Cosh »L/2 =1)) (8) 
Cosh AL/2 


Ww, =transverse frequency with coupled longitudinal motion 
=transverse frequency as calculated disregarding the offset 


w 
oe of the center of gravity from the center of stiffness 


te =total deflection of the assumed mode shape 
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Y =first approximation to the additional deflection 


V5 = second approximation to the additional deflection 


When substitutions shown in Appendix A are made, all the terms in brackets 


in equations (7) and (8) are, to an approximation, dependent only on fre- 


quency ratio, f. 


Yot 
f= a: 
of 
where 
iA He longitudinal frequency as calculated disregarding the 
offset of the center of gravity from the center of stiffness. Specifically, 
2 
W 2 
2 of 4.73 
pa» = () YX) 
ikae Wy © Lon 
pa? = (4:13) ALY 
1 ~ NB in 
Similarly, 


L/2 _ + 


1 


rja 


Another important result derived in Appendix A is that the frequency ratio, 
f, is independent of mass per unit length. 

These two important results suggest a rearrangement of the 
frequency equation 4, to permit easy calculation of the change in frequency 


for any system. 


e 
it) .— 
“ot “ae Wee 
3 XY 


> = 
= | = — [ae 
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This equation can be put in a still more usable form. 





et y x (9) 
Wo =2 =e 2 9 
Note AL g(r) (AL) ale) 


where g(f) and h(f), functions which are only dependent on frequency ratio, 


can be obtained from Figures I and II. 


The function for the first approximation to beam deflection is: 


) 
32. ee) (Cos #3 - 1) 


Cos AL Lewin Mame. 32 Cog aie - 1) , (Cos ne - 1) 1 
wee Ge =) )cos (£5) Kc -(42)° ) Cos os (a + Ass =)*) (cosh 21/2) 


[2 Cos (AL/2 - 1 - 





The function for the second —" to beam deflection is: 


h 
(f) ea - ats <a Me oy 8 2). @ ) 
h(f) = Cos AaL/2 oe on 
; 1- (ey) (a - Gy") cos 1 - 


TT 
ve (20) ee P , {Cos »AL/2 = 1) 
2 


(Cos £n)* (1 - (22)*)* Cos AL/2 


-~ va 
, (Cosh »L/2 - 1 ) ]=[e See NT Bo ak 


aS (A2)*)° Cosh ee 


Cosh AL/2 
These expressions do not have to be evaluated in finding the change in 
natural frequency. All that has to be evaluated is Wot? wo and = with 
g(f) and h(f) from Figures I and II to use the approximate analytic method. 
Above a frequency ratio of one-half the approximate analytic method 
gives poor results. However, it is still desirable to estimate frequency 
changes at higher frequency ratios. An examination of the relative order 


of magnitude of these two functions suggests an alternate method of determining 


frequency changes from data of only one numerical example of the lumped mass 
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FIGURE I 
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FUNCTION FOR TSE SECOND APPROXIMATION 
TO THE CHANGE IN BEAM DEFLECTION 
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method shown in Figure IV. This alternate method is based on the fact that 
2 
h(f) ~ a(t) | 


To use Figure IV for other beams use the following procedure: 


W 
AG Determine — from this figure for the desired frequency ratio. 
ot 


2. This number can be used to solve Equation (10) for x. The 
solution is shown as Equation (11). 


=2 
3. Use this value of x in a new frequency equation with 


for the desired application. See Equation (12). 


(10) 


(12) 


(12) 





where 
“e 
oe frequency change obtained from Figure IV 
ot 
=2 
- = . 
77 wr 
Ay= 
a = for the desired application 
w.° 
— = frequency change for the desired application. 
ot 


B. Results of the Lumped Mass Method. 


The equations for determining the conditions at each mass along 
the beam are shown below. The zero subscript refers to the conditions at the 
previous section; the 1 subscript refers to the mass at the point being 


evaluated. The actual numerical value of each of these equations becomes the 
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subscript zero term when the next point along the beam is evaluated. 


s = 2 = 
M>=M +Ve+y m, Ww (u, - y 8, ) 


2 
vy 3 “ “we Fy 





Ms vis 
1 °o Et * Sen 
Me2- vig 
al = Vo i‘ we i a hi eran 
re 
mer ak: Ke 


2 a 
FL oe F,+m,w (u, - ¥ 6) 


where 


bending moment 


shear force 


i 


longitudinal force 


M 
V 
F 
K 


spring constant between masses 


fa distance between masses 


One advantage of the lumped mass method is that it is useful in deriving a 
frequency equation. For the three mass systems shown in Figure III, the 


frequency equation is: 


~2 yt -2 3 
y iL 
O+2-mw (5 ir + oes 2) + mew! (toe + 3 om) 


C. Results of a Numerical Example. 


In order to check the validity of the results, a numerical example 
was done using both methods. The change in frequency is shown plotted in 
Figure IV as a function of frequency ratio, f. The plot is for a system of 


these properties: 
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A= 1 in 
mM 

t= sae in 

y = 1/h in 

a = .284 lbs/in? 


=2 
ak? 


In the lumped mass method the system was divided into eleven masses (Figure III). 
For system of fewer masses the agreement between methods was poor. For example, 


at a frequency ratio of one-fourth: 


W 

— = .97h analytic method 
ot 

"t 

— = -961 eleven mass system 
ot 

“4 

a = -950 nine mass system 
ot 

ue 

aT -—910 seven mass system 
ot 


For all numerical work, Be and Woy Were computed from: 


22.4 if 
Wot ~ = \/ any 
_ EgA 
- V (aA) 


D. Numerical values for Marine Application. 
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Vos 





The approximate analytic method indicates that there are two sig- 


nificant parameters in calculating the change in transverse frequency, namely, 


=2 
frequency ratio, f, and the geometric factor — . The geometric factor 


has been evaluated for two ships: 





#16= 


2 
merchant ship (3) — = 087 


submarine (4) a ~ .O17 


These numbers, of course, are dependent on the condition of loading and, 
furthermore, they are valid only for the midship section. The geometric factor 
for the merchant represents an upper limit of this number for the "Gopher 
Mariner". Approximate frequency ratios for these ship types are: 

submarine (4) fx .15 


merchant ship fx .25 





ie 
DISCUSSION OF RESULTS, CONCLUSIONS, AND RECOMMENDATIONS 


A. Accuracy of Methods 
The approximate analytic method is inaccurate above the frequency 


ratio of one-half. The lumped mass method, on the other hand, can be used 
over the entire range of frequency ratio. 

Reference (5) states that to ensure better than one percent 
accuracy in computing frequencies by the lumped mass method, thirteen 
masses should be used per wave length. If more masses had been used in the 
numerical example, the agreement between methods at the low frequency ratio 
would be better than it already is. For eleven masses there is only a 1.3% 


adifference. 


B. Significance in Ship Vibrations. 
When the constants for the submarine and merchant ship stated in 


Part III above are used in Equation (9) in conjunction with Figures I and II, 
the change in transverse frequency is truly insignificant for the fundamental mode. 
For higher modes of vibration, the frequency ratio is larger and the change 


in frequency might not be so insignificant. 


C. General 

Two general methods for calculating frequency change due to an 
offset of the center of gravity from the center of stiffness are given . 
One method is approximate; the other is exact. Also, a method of estimating 
frequency change by a simple calculation from one lumped mass numerical 
example is explained. The coupling of the transverse and longitudinal modes 
of vibration in any structure can be estimated or obtained exactly by the 


methods explained in the text of this thesis. 





Soe 


All the work is restricted to the fundamental mode of vibration. 
The lumped mass method, of course, could be used for higher modes. It is 
recommended that the coupling due to the offset of the center of gravity 


from the center of stiffness for higher modes be investigated. 
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APPENDIX 
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APPENDIX A 


DETAILS OF THE APPROXIMATE ANALYTIC METHOD 


The approximate analytic method, outlined below, is divided into 
steps, the numbering of which corresponds to the numbering in the procedure 
Section IIA above. The description which appears in that section is not 
repeated. 

Step 1; Derivation of the Equations of Motion 
With the sign convention shown in Figure A-l an equilibrium 


force balance on a particle of the system yields this differential equation 


of motion. 
2 ys 
- SoA) Se (u+ Zz) + Org -0 (1) 
at ax 
where: 


(aA) = weight per unit length 


gravity 


08 
i 


y= a= s, = distance of center of stiffness from the center of 


gravity 
E = Young's modulus 
u = displacement in the x- direction of the center of stiffness 
from the equilibrium position 
du 
S=- ES 
2 
: 2 
Letting - on. C= wet) 
ot 
the equation can be rewritten 
2 
Aloe ( a cee 7) + EACLE = 0 (la) 
g dg ax 


A moment balance about the center of gravity of & small particle 


can be written using the sign convention of Figure A-II. 





Jone 
FIGURE Al 


LONGITUDINAL FORCE LALANCE 
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Yo > distame to the center or bending 


ee = Gistance to the center of gravity 


ds; 
(3 + a5 ax )A 
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FIGURES A-II 
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V ax+M - (M + am) + (9 +S ax - 8) ay = 0 (2) 


Let . 
av (ef) SY ax 


BE — 


ca 
" 


A = area of the cross section 


Equation (2) can be rewritten: 


2 4 3 
QA oy , prof + ma SUF - 
© st" ox sx2 


Step 2: Assumed Mode Shape wwii! Mourns 
According to the classical theory olf vibraticn, the equation of 


transverse motion for a uniform beam free in space is: 


Cos AL/2 
y = e| Cos AX + ane 5) Cosh Ax_| (3) 
L = length of the beam 


x = distance from mid-length 


eo) 
i 


arbitrary constant 


For the fundamental twoenoded mode of vibration: 


see 
AL = 4.73 = Vv, we a L (3a) 


Step 3: Longitudinal Motion 
Equation 3 is substituted into Equation 2 and the resulting 


differential equation solved for the longitudinal motion, u. The arbitrary 
constants are evaluated by assuming no longitudinal motion at mid-length 


and no stress at the ends. 





api. 


xa= 0 Oe 
du 
x =L/2 ae! 


The equation for longitudinal motion becomes: 


= F Sin ~—+ P_ Sin »4x + F. Sinch rx (4) 


aL VD, 2 S 


where 
2 D, »* 
Fo= - D os AL/2 
1 AE Ds ines au iv joy | 
D 
aa : D 
1 «= DS Xr 
D,~ D 
we. ae 3 
3 1+D Ne 
_ SBA ( 
a 5) 
1 we (aA) 
D,=ye 2d (6) 


- Cos AL/e 
Dy ° - Sabie a 


Step 4: Additional Deflection Caused by Longitudinal Forces 





In this step the objective is to find the additional vertical 
deflection caused by the longitudinal motion and forces. Again referring 
to Figure A-2, the differential equation to be satisfied is: 


OM = 
ae * a Ay = 0 (8) 


where 


additional transverse deflection 


md 
it 





a 


Equation (8) becomes 


yt a-s- au 
3 es 7 Yaa (9) 
ax ax 


These boundary conditions are applicable in the solution of equation (9): 


a y' 

— = re) yosab/2 
dx 

du 

== =e x = L/2 
g 

wv = 0 <= 
u=0QO x = 0 


The result of substituting Equation (4) into Equation (9) and 
solving the differential equation is the equation for the additional deflection 
of the beam caused by the coupled longitudinal motion. The equation is: 


F 


yi = [Fy yD, < ~ 2 ] 
=> |F Cos i + Cos Ax - 5 Cosh Ax} + Constant (10) 


A measure of the total additional deflection in the beam is obtained by 


evaluating Equation (10) at midlength and at the end. 


- vy! es g 
Y1 = 9 x20 Vy = 1/2 (11) 
202 4h (Cos pot 2) 

_ Ay \ (Cos a /2 = 1) 

Y, == -¢Cos xL/2 oe wae 

" Cos = (1 = D, 2) Cos AL/2 
Pa VD, 

Cosh AL/2 = 1 

es (12) 
1+D, ,2 Cosh aL/2 


Step 9: The Second Approximation to Total Transverse Deflection: 





When Equation (10) is then cycled through Steps 3 and 4 and 
evaluated at the mid-length and end point, the result is a second approximation 


to additional deflection of the beam. 





whe 


| L/2 
cee 2D,” a" a -, an 
= (=) (y) ¢ Cos AL/2 ae a nes a L/2 
D 
1 


ab ee 


=f ) Sin 1/2rPa_ (Cos AL/2 - 1 
-_) (He L al Cos AL/2 
1-D oe 
(cosh XL/2 - 1 
een: gosh ELE = 2) (13) 


Step 6: Calculation of Natural Fre 





With an assumed frequency and deflection, the Stodola procedure 
is to calculate the inertia loading and resultant deflection. If the 
assumed frequency equaled the natural frequency, the resultant deflection 


would equal the assumed deflection. Mathematically: 


wy Y= Wo Jo (14) 
_ WS oo (15) 
Sear assumed deflection 

yi = resultant deflection 

et ae assumed frequency 

wae natural frequency of the beam 


When the Stodola equation is applied to the approximate analytical 
procedure, the equation for determining the natural frequency of a system 
with a coupling of longitudinal and transverse motion becomes: 


a 


2 2 re) 
t ot a + Yon t Y,, 


Yo is Equation (3) evaluated at midlength and at the end 


Y, is Equation (12) 


Y, is Equation (13) 





soy 


x 
It 


a natural transverse frequency of the system without coupling 


of longitudinal and transverse motion. 


x 
tt 


‘ natural frequency of the system with coupling of longitudinal 


and transverse motion. 


Additional information can be obtained by some algebraic 
manipulation. From Reference (2), the equations for the square of the 
natural longitudinal frequency and transverse frequency of a uniform beam 


are respectively: 





(w)2 = @) Tay 
of L aA 
4° 
(w IF - Ce. Egl 
ot L (oA) 
Bierefore, a 
fe = Not _ 22..4o aL (16) 
7 Qo 2 2 A 
Woe TT L 


Equation (16) gives the very important result that the frequency ratio, f, 
is independent of mass per unit length. 


An examination of the expression D he yields another important 


iL 


result. Equations (3a) and (5) show that D, Ne Ae equal to: 


D, 0° = ts oA | (+13 13) 





18 
ai 


ro ae 


2 
W 
Lee > 
(24) (8) 
Y% 


If the change in transverse frequency is small, Wie WO 





Sage 


Therefore, 


ut 


pr»? = 3) (4B) 


Similarly, it can be shown that: 


L/2 a 
i Se 


These manipulations show that Equations (12) and (13) are made up to two 
-2 


factors; one is -f and the other is a complicated expression which is 


dependent only on frequency ratio. 
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APPENDIX B 


DETAILS OF THE LUMPED MASS SYSTEM 


The equations in the beam sequence calculation are derived from 


the fundamental equation for bending in a bean. 


In the case of a beam with the center of gravity and center of stiffmess 
coincident, this equation and its integrals, using the sign convention 


of Figure B-I, become: 
2] 





BI 3 =m +V x 
Abe O O 
2 
cy x 
ee ax aon sf Ve 2 Bacal 
Mx" Vix 
Ely =-3 nee ca SS 
where 
ve = shear at thc boundary or beam end 
M, = moment at the boundary or beam end 


When these quations are used in a step-by-step procedure, the equations 
for moment, shear, slope, and deflection when evaluated at an incremental 


distance to the right of mass 1, shown in Figure B-lI, become; 





M=M +VE 
Ve=V + wey 
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FIGURE B-I 


SIGN CONVENTION FOR TNE LUMPED MASS SYSTEM 
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where 


deflection at the end 


ed 
it 


@ = slope at the end 


In the case of where the center of gravity and center of stiffness 


are not coincident, these equations are modified and two additional equations 


are added. 


a a z 
M = M. + Ve +ym Ww (u, 2 ¥y 8, ) (2) 
V. =V_ + Ww y (2) 
i eee 1 
Mt Vie 
Gin o. @EL SET (3) 
Me” ve 
¥y = Yq * Sot * oer * CET “) 
F 
o. os 
* (uy - 5S (6) 
ae PF +m, Ww (u, ~ ¥ ey) 
where 
Ce longitudinal displacement at position zero 
u) = longitudinal displacement at position one 
K = spring constant = 
ye = longitudinal force at the end of the beam 
F, = longitudinal force at position 1 
y = distence of the center of gravity from the center of stiffness 


In the case of a beam vibrating free in space, the conditions at the end 


of the beam are further specified. 
2 z 
= ae (u, - y 6.) 


” D7. = 
M,=FLy =m, y (uy -y 8) 


< 
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2 
mW Y, 
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Equations (1) through (6) can be evaluated at other positions along the beam 


using a table as shown in Appendix C. To start the calculation, let 


eo 
ue=uUu 
O Oo 
9 = @ 
oO oO 


Then a table is completed for each mass along the system using the results 
of the previous table to start the next. The end of the last table, which 
will be the conditions at the mid-point of the beam, will yield three 
equations: 


9 = ayy ~ ao) Uo a7 seal 80 


B15 + B55 Ww, + B25 oe 


Uu= 83 + 253 Uy + a2 hs 
If the proper frequency was used in the calculations, these equations all 
equal zero. The procedure is to let u and V equal zero, solve for Us and 


86 and substitute into the @ - equation. When § equals zero, the proper 


frequency was used in the calculations. 
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APPENDIX C 
CALCULATIONS OF LUMPED MASS METHOD AT A FREQUENCY RATIO 


OF .25 WITH ELEVEN MASSES 


ae 4) in” 
1 = 1/12 in.” 
a = .264 Ibs./in.? 
L = 8.2h in. 
f «= 25 
y =1/4 in. 
L = .02h in. 
m = .60666 x 10° (1bs.-sec.“/in. ) 
Woe = 3.69543 x 10° (radiens/sec. )* 
at = .03295 x 10°? (1/lbs.-in) 
ar. .01357 x 107 (i/1bs.) 
PEI 
4 “5 
Sez = -00372 x 10 (in. /lbs.) 
2 = .O00274 x 107? (in 1bé. ) 
assume w= .975 We 
Therefore, 
mw° = 2.13120 x 10° 
F, = 1.06560 x 10? u, = 26640 x 10 9. 
M, = .26640 x 10? u, ~ .06660 x 10° 6, 
v5 = 1.0650 x 10° y. 


The five-page 


the following 


tabular calculation to reach the center of the beam is on 


pages. 
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The last page of the tabular calculation yields three equations which would 


all equal zero if the right frequency was assumed at the beginning of the 


calculation. 
8 = 1.27136 yee 23484 aes 2.14740 0, 
V = 13.68628 aa -88648 aes 23.55147 Q, 
u = .00219 + AS + .92853 Us + .O1907 on 

let ie i 
ane) 
V=0 

Solve for U, and os and substitute into the §@ equation. 
Uo : -00958 
0, == 58148 
@ = .02493 


When the assumed frequency is .95 Wot? and the five-page tabular calculation 
repeated, 

@=- .01893 
The result is plotted in Figure C-l,and thus the natural transverse frequency 
with offset center of gravity is: 


W= 961 oe 
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